First-principles calculations of magnetization relaxation in pure Fe, Co, and Ni with 

frozen thermal lattice disorder 
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The effect of tlie electron-pfionon interaction on magnetization relaxation is studied within the 
framework of first-principles scattering theory for Fe, Co, and Ni by displacing atoms in the scatter- 
ing region randomly with a thermal distribution. This "frozen thermal lattice disorder" approach 
reproduces the non-monotonic damping behaviour observed in ferromagnetic resonance measure- 
ments and yields reasonable quantitative agreement between calculated and experimental values. It 
can be readily applied to alloys and easily extended by determining the atomic displacements from 
ab initio phonon spectra. 

PACS numbers: 72.25.Ba, 72.25. Rb, 72.10.Di 



I. INTRODUCTION 

The drive to increase magnetic storage densities and re- 
duce access times is focussing renewed attention on mag- 
netization dynamics in response to currents and external 
fields.^ The dynamics of a magnetization M in an effec- 
tive field Heff is usually described with the phenomeno- 
logical Landau-Lifshitz-Gilbert equation 



dM 

dt 



= -7M X Hcfi' + M X 



A(M) rfM' 
7MI IT 



(1) 



where Ms = |M| is the saturation magnetization den- 
sity and 7 — g^iofJ-B/fi- the gyromagnctic ratio expressed 
in terms of the Lande g factor and the Bohr magneton 
/iB- The first term describes the precessional motion of 
the magnetization in the effective field that includes the 
external applied field, the exchange field, anisotropy and 
demagnetization fields. The second term describes the 
time decay of the magnetization precession, the Gilbert 
damping,^ in terms of A(M) that is in general a symmet- 
ric 3 x 3 tensor.'^ For isotropic media, the damping is fre- 
quently expressed in terms of the dimensionless parame- 
ter a given by the diagonal element of A, a = X/jMg. 

There is general agreement that spin-orbit coupling 
(SOC) and disorder are essential ingredients in any de- 
scription of how spin excitations relax to the ground 
state. In the absence of any other form of disorder, 
one might expect the damping to increase monotonically 
with temperature in clean magnetic materials and in- 
deed, this is what is observed for Fe in ferromagnetic 
resonance (FMR) measurements.''^'^ Heinrich et al.^ de- 
veloped an explicit model for this high-temperature be- 
haviour in which itinerant s electrons scatter from local- 
ized d moments and transfer spin angular momentum to 
the lattice via SOC. This s-d model results in a damp- 
ing that is inversely proportional to the electronic relax- 
ation time, a ~ l/r, i.e., is resistivity-like. However, 
at low temperatures, both Co and Ni exhibit a sharp 
rise in damping as the temperature decreases. The so- 
called breathing Fermi surface model was proposed*^^'^ 



to describe this low-temperature conductivity-like damp- 
ing, a T. In this model the electronic population lags 
behind the instantaneous equilibrium distribution due to 
the precessing magnetization and requires dissipation of 
energy and angular momentum to bring the system back 
to equilibrium. 

Of the numerous microscopic models that have been 
proposed^^ to explain the damping behaviour of metals, 
only the so-called "torque correlation model" (TCM)^^ is 
qualitatively successful in explaining the non-monotonic 
behaviour observed for hep Co. An effective field ap- 
proach can be used^'^^^^ to identify conductivity-like and 
resistivity-like behaviour at low and high temperatures, 
respectively with intraband and interband terms in the 
TCM." Evaluation of this model for Fe, Co and Ni us- 
ing first-principles calculations including SOC for the 
host electronic structure and a band-, wavevector- and 
spin-independent relaxation time approximation (RTA) 
to model disorder yields results for the damping a in good 
qualitative and reasonable quantitative agreement with 
the experimental observations. ^^"^^ The disadvantage of 
the RTA is that it is difficult to unambiguously map mi- 
croscopically measured disorder onto a unique value of 
the relaxation time r. 

A formulation of magnetization damping in terms of 
scattering theory, that is equivalent in linear response to 
the Kubo formalism, was recently applied to the study 
of substitutional alloys with intrinsic disorder yielding 
good agreement with experiment without introducing 
any parameters."'^'' The discrepancies remaining between 
the experimental data measured at room temperature 
(and higher) and the T = calculations pose questions 
about the role of various types of thermal disorder. In 
this paper we combine the scattering theory formulation 
of Gilbert damping with structural lattice disorder to 
model finite temperature lattice effects. Our main result 
is to show that this can reproduce the non-monotonic be- 
haviour of the magnetization relaxation as a function of 
temperature. 

The paper is organized as follows. In Sec. II, the 
"frozen thermal lattice disorder" scheme is introduced 
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and a brief description of the scattering theory is given. 
This is followed by some technical details of how the 
calculations are performed and how the resistivity and 
Gilbert damping parameter are determined. The results 
are presented and discussed in Sec. Ill and compared to 
previous calculations and experiments. A short summary 
and some concluding remarks can be found in Sec. IV. 
In the Appendix, we discuss a factor of An commonly 
omitted when Gilbert damping frequencies are given in 
Gaussian units. 



II. THEORETICAL METHODS AND 
COMPUTATIONAL DETAILS 




FIG. 1. (Color online) Schematic picture of the proposed 
frozen thermal lattice disorder. Atoms (red dots) on the ideal 
lattice (left panel) are displaced with a random Gaussian dis- 
tribution and form a static configuration (right panel) in the 
electronic transport calculation. The displacement of atom i 
is denoted as Ui (blow up). 



Assuming the Born-Oppenheimer approximation,-'^® 
static disorder is introduced in the transport calculations 
by displacing atoms rigidly and randomly from their ideal 
lattice positions in what we call a "frozen thermal lat- 
tice disorder" scheme. In the temperature range we are 
interested in, far below the melting point, typical dis- 
placements are of the order of several hundredths of an 
angstrom, small compared to the lattice constant. We 
can therefore adopt a harmonic approximation and cor- 
responding Gaussian distribution of displacements char- 
acterized by the root-mean-square (RMS) displacement, 
A = ( I Ui p ) where the angular brackets indicate an av- 
erage in which the index i runs over all atoms (Fig. 1). As 
the temperature increases, higher energy phonon modes 
are occupied so A increases. In the present study, this 
qualitative correlation between temperature and A is suf- 
ficient to produce a non-monotonic damping. If we knew 
the phonon dispersion relation from a Debye model or 
ab initio calculations, the atomic displacements could be 
determined explicitly by summing contributions from all 
vibrational modes occupied at a specified temperature. 
Such a description of the lattice disorder introduced by 
finite temperatures could then be straightforwardly com- 
bined with scattering theory to study temperature de- 
pendent magnetization relaxation quantitatively. 

It was shown by Brataas et al}^ that, for a single do- 
main ferromagnetic metal (FM) sandwiched between left- 
and right-hand leads of non-magnetic (NM) material, the 
Gilbert damping tensor G can be expressed as 



Gij(m) = A(m) 



V = ^Re 



Tr 



dS dS^ 



dm; drrii 



, (2) 



where V is the volume of the ferromagnet, the scattering 

f r t' \ . . . 
matrix = ( ^ r' ) '^^ given in terms of reflection and 

transmission matrices for Bloch waves incident from the 
left (r and t) or right (r' and t') leads. When SOC is 
included, S depends on the direction of the magnetiza- 
tion unit vector m = M/Afg. The microscopic picture 
of magnetization damping implicit in the scattering for- 
mulation is of energy being transferred slowly from the 
spin degrees of freedom through disorder scattering and 



SOC to the electronic orbital degrees of freedom and then 
being rapidly lost to phonon degrees of freedom in ther- 
mal reservoirs attached to the leads. From the transmis- 
sion matrices, we can also calculate the conductance of 
the system within the Landauer-Biittiker formulation as 
G= (eV/i)Tr{ttt}. 

To evaluate the scattering matrix at the Fermi level, 
we use a "wave-function matching" scheme^^ imple- 
mented with tight-binding linearized muffin-tin orbitals 
(TB-LMT0s)2° that was recently extended to include 
SOC." The electronic structure of the NM|FM|NM sand- 
wich is first determined self-consistently using a surface 
Green's function method^^ with a minimal basis of TB- 
LMTOs in the atomic sphere approximation. In the 
current study of Fe, Co and Ni, we consider Au|Fe|Au, 
Cu|Co|Cu and Cu|Ni|Cu sandwiches so the lattice con- 
stants of the NM leads and FM scattering regions match 
almost perfectly.^^ The two-dimensional (2D) Brillouin 
zone (BZ) of the 1x1 unit cell is sampled with a 120 x 120 
grid in the self-consistent calculations. 

Disorder is introduced by randomly displacing atomic 
spheres in the FM scattering region using the frozen ther- 
mal lattice disorder scheme described above. The calcula- 
tions are rendered tractable by imposing periodic bound- 
ary conditions transverse to the transport direction. It 
turns out that good results can be achieved even when 
these so-called "lateral supercells" are quite modest in 
size. In practice, a 4 x 4 lateral supercell and a 28 x 28 
2D BZ grid were found to be sufficient for Fe, and a 5 x 5 
supercell and a 32 x 32 grid for Co and Ni, respectively.^'^ 
The thickness of the FM region ranges from 20 to 340 
atomic monolayers. For every thickness of the ferromag- 
net, we average over a number of random disorder config- 
urations. The sample-to-sample spread is small for large 
values of A and five configurations are sufficient; for small 
values of A, as many as 35 configurations are used. 

For a ferromagnetic slab of thickness L, we write the 
resistance of the system as R{L) = 1/G{L) = 1/Gsh + 
2i?if -|- Rh{L), where G{L) is the total conductance, and 
Gsh the Sharvin conductance of the ideal leads; 
is the NM|FM interface resistance and Rh{L) the bulk 
contribution. ^^'^^ When the ferromagnetic slab is suffi- 
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FIG. 2. (Color online) Total resistance (a) and damping (b) 
as functions of the thickness of the ferromagnetic slab for 
Au|Fe|Au(001) with A = 0.0259 apc. Dots indicate the calcu- 
lated values averaged over five configurations while the solid 
line is a linear fit. 



ciently thick, we expect to recover ohmic behaviour with 
Rh{L) w pL; this was demonstrated explicitly for the 
case of alloy disorder in Ref. 17. As shown in Fig. 2(a), 
this expectation is borne out by the present calculations 
and the resistivity p arising from the frozen thermal lat- 
tice disorder can be extracted from linear fitting. We 
write the damping parameter G analogously as a sum 
of an interface contribution Gk and a bulk contribu- 
tion, Gh{L). If we further express the bulk contribu- 
tion in terms of the dimensionlcss damping parameter 
a as Gh{L) = X ■ V = ajMgAL where A is the cross 
section, then we expect the calculated damping to grow 
linearly with the thickness of the ferromagnetic layers, 
G{L) = Gif -I- a'^MgAL. Once again, this expecta- 
tion is borne out by the calculations as demonstrated 
in Fig. 2(b). and a can be determined from the slope. 



III. RESULTS AND DISCUSSION 

The resistivities and damping parameters calculated 
for bulk Fe, Co and Ni are shown as a function of the 
RMS displacement in Fig. 3. For all three materials, p 
increases monotonically with A as expected and a agrees 
with the predictions of the torque correlation model. ^'^"^^ 
For large values of A, a calculated for Fe and Co is found 
to increase with increasing A, while tending to saturate 
for Ni in agreement with experiment. ^'^'^ For small values 
of A, a increases rapidly as A decreases. This sharp 
rise is the conductivity-like behaviour observed at low 
temperature for Co and Ni, demonstrating that a simple 
model of frozen thermal lattice disorder can reproduce 
the non-monotonic Gilbert damping seen in experiment. 

To the best of our knowledge, this work represents the 
first calculation of p for these ferromagnetic metals, only 
nonmagnetic materials having been considered in pre- 
vious first-principles studies of temperature dependent 




FIG. 3. (Color online) Calculated Gilbert damping and resis- 
tivity for bcc Fe, hep Co and fee Ni as functions of the RMS 
displacements measured in units of the corresponding lattice 
constants, a. The error bars reflect the configuration spread. 
Experimental resistivities^* are used to label a number of re- 
sistivity values with a temperature. 



resistivity.^^ While there is ordcr-of-magnitudc agree- 
ment with experiment, we cannot make a rigorous com- 
parison if we only have a qualitative knowledge of how 
A depends on temperature T. If we adopt the com- 
monly made assumption that the effect of temperature 
on p and a can be expressed in terms of a phcnomcno- 
logical scattering time r(T), i.e. p{T) = p{t{T)) and 
a{T) = a(r(T)), then a should have a well-defined de- 
pendence on p. In the present case, this amounts to 
assuming that our lattice disorder can be mapped onto 
t(T) = r(A(r)). We can make a first qualitative com- 
parison by using the experimental p{T)'^^ to indicate a 
number of temperatures in Fig. 3. Around these temper- 
atures, the corresponding damping parameters behave in 
the same way as observed in expcrimcnt:^'^''' the damp- 
ing for Fe has a minimum around 273 K and increases 
slightly at about 500 K; at 273 K, the damping for Co has 
increased away from its minimum value; while that for Ni 
has already saturated. This qualitative agreement indi- 
cates that the correlation between p and a is reasonably 
parameterized by A. In Table I wc compare the mini- 
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TABLE I. Minimum values of the Gilbert damping A in units 
of 10** s^^ with respect to: temperature (experiments); re- 
laxation time, T (Torque-correlation model); and A (present 
work). Experimental damping frequencies have been multi- 
plied by 47r; see the Appendix for details. 



A 


bcc Fe [001] 


hep Co [0001] 


fee Ni [111] 


Expt. 


8.8'', 4.8'" 


ga 


29", 28" 


TCM^^ 


5.4 


3.7 


21 


This work 


3.9 


2.3 


20 



"Ref. 5; ''Ref. 4; "Ref. 7. 



mum values of damping calculated as a function of A with 
the corresponding minimum values calculated vi/ithin the 
TCM as a function of the relaxation time r reported in 
Ref. 13 and also with the minimum values determined 
experimentally as a function of temperature.'*'^'^ We see 
that our minimum values are lower than the TCM values 
that are, in turn lower than the experimental values. 

There are two noteworthy qualitative differences be- 
tween our results for small values of A and the low- 
temperature experimental observations, (i) For Fe, we 
find a conductivity-like damping behaviour that has not 
been seen in experiment^'^ but has been found in the 
TCM calculations. ""^^ (ii) For Ni and Co, the satura- 
tion of the damping observed at very low temperatures^''' 
is found neither in our calculations nor in the TCM 
studies. ^'^"^^ We suggest that these discrepancies arc in 
part related to extrinsic contributions to the Gilbert 
damping that compete with the intrinsic thermal effect. 
There will always be some amount of extrinsic disorder 
such as impurities and defects in experimental samples 
whose contribution to damping is essentially tempera- 
ture independent and becomes dominant when the ther- 
mal disorder becomes negligible at sufficiently low tem- 
peratures. If the extrinsic disorder in the Fe samples 
measured in Refs. 4 and 5 was so high that Fe was in the 
resistivity-like regime^'' for all temperatures, then reduc- 
ing the extrinsic disorder would be expected to lower the 
damping. This is consistent with the minimum damping 
value we calculate being smaller than the experimental 
value, see Table I, since our calculated damping comes 
only from the electron-lattice scattering, without any ex- 
trinsic contributions taken into account. For Ni and Co 
that exhibit both conductivity-like and resistivity-like be- 
haviour, the discrepancies cannot be explained away so 
simply. Better characterisation of the experimental sam- 
ples as well as further theoretical study are required. 

Comparing the minimum values of A, as we have just 
done, assumes that t (or A) characterizes A and p com- 
pletely. In particular, it assumes that the effects of var- 
ious temperature-dependent disorder scattering mecha- 
nisms can be represented by a simple relaxation time. 
If this is correct, then the microscopic details of the 
disorder that give rise to particular values of p and A 
should not be of paramount importance and we should be 
able to use experimentally determined {X{T),p{T)} and 
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FIG. 4. (Color online) Gilbert damping frequency as a func- 
tion of resistivity for bcc Fe, hep Co and fee Ni. Calculated 
results are shown as lines: for frozen thermal lattice disorder 
as (black) solid lines and for frozen spin disorder as a (red) 
dashed line (Fe only), g = 2 is used in the calculations. Sym- 
bols are experimental damping values from Ref. 5 (0)i R-^f- 4 
(□), and Ref. 7 (A). The corresponding resistivity values are 
take from References 25 and 28. Experimental damping fre- 
quencies have been multiplied by 4n; see the Appendix for 
details. 



theoretically determined {X{A) , p{A)} values to compare 
•^oxpt(p) and Acaic(p)- Because we are not aware of any 
simultaneous measurements of A and p, we have com- 
bined A(T) and p{T) from different experiments*'^''''^^'^^ 
to plot A(p) in Fig. 4; the calculated results are shown as 
black solid lines while the available experimental data are 
shown as symbols. To facilitate the comparison, the cal- 
culated damping parameter a in Fig. 3 is converted (see 
the Appendix) to damping frequency A using A = a'^M^. 
For convenience we have used g = 2 for the calculations 
instead of the measured ^-factors. 

The best agreement between the calculated curve and 
Bhagat and Lubitz's experimental values^ is obtained for 
Ni. As mentioned above, including the scattering associ- 
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ated with residual resistivity should reduce the damping 
calculated at low temperatures (resistivities) and increase 
it at high temperatures (resistivities) giving better agree- 
ment with experiment. However, it is not clear how the 
later measurements by Heinrich et al/ might be accom- 
modated in this picture. Clearly, there is a need to deter- 
mined which experiment most accurately represents the 
"intrinsic" damping in Ni. A similar situation obtains 
for Fe where there is an even larger discrepancy between 
the two existing experiments. The figure highlights the 
disagreement between a conductivity-like behaviour at 
low values of resistivity in the calculations and its ab- 
sence in both measurements. What is also striking is the 
failure of the calculations to reproduce the high temper- 
ature (resistivity) enhancement. For Co, where there is 
only one measurement, we find that in a reversal of the 
situation for Fc the conductivity-like behaviour is much 
more pronounced in experiment than in our calculations. 
It is unclear from this and the TCM studies how much of 
these discrepancies might result from subtle features of 
the electronic structure that are described inaccurately 
by the local spin density approximation. Tests with dif- 
ferent exchange correlation potentials and increased max- 
imum angular momentum cutoff indicate that the corre- 
sponding changes to the electronic structure at the Fermi 
energy cannot explain the discrepancies. 



We can test the uniqueness of t directly by calculating 
\{p) for different types of disorder that give rise to the 
same resistivity. We do this for Fe by modelling frozen 
thermal spin disorder in a manner analogous to the way 
we have modelled frozen thermal lattice disorder. We 
introduce a random, Gaussian distribution of spins with 
respect to 9, the polar angle between the local atomic 
magnetization and the global magnetization directions, 
together with a uniform random distribution in the az- 
imuthal angle 0. Together, 9 and (j) determine the local 
magnetization direction. Wc then calculate A and p as 
a function of the root mean square polar rotation an- 
gle. The results are included in Fig. 4a as a dashed (red) 
line. Though the qualitative behaviour of Aiattice(p) and 
Aspin(p) is the same, quantitatively they are clearly dif- 
ferent, differing by as much as a factor of two for the 
lowest value of resistivity shown. This is clear direct ev- 
idence that different microscopic scattering mechanisms 
contribute to resistivity and Gilbert damping differently 
and that it may not be sufficient to use a single electronic 
scattering time to represent various types of disorder. It 
is also of interest to study to what extent the effect of dif- 
ferent scattering mechanisms are additive. To do so, it 
is desirable to introduce "real" temperatures by calculat- 
ing the phonon and magnon spectra from first-principles. 
This is a major undertaking and will be the subject of a 
separate study. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we report the results of first-principles 
calculations of the Gilbert damping with thermal lattice 
disorder for pure Fe, Co and Ni in the framework of a 
recently introduced scattering theory. The effect of tem- 
perature on the lattice is simulated by displacing atoms 
with a random, Gaussian distribution. Our main result is 
that both the conductivity- like and resistivity-like damp- 
ing behaviour observed in FMR measurements are repro- 
duced by the scattering theory. The reasonable quan- 
titative agreement between our results and experiment 
demonstrates that our simple thermal disorder scheme 
accounts for the dominant, intrinsic effect of lattice tem- 
perature in magnetization relaxation. By calculating the 
damping as a function of resistivity and comparing the 
results to experiment, we highlight discrepancies between 
different experiments and between the calculations and 
experiments. Part, but not all, of the discrepancies for 
the Gilbert damping can be attributed to competition be- 
tween extrinsic and intrinsic scattering. An exploratory 
calculation of thermal spin disorder for Fe indicates that 
different types of disorder affect the Gilbert damping and 
resistivity in different manners. This work needs to be 
extended to different materials before more general con- 
clusions can be drawn. It is of particular importance to 
study how different types of disorder combine to affect 
the damping. To avoid having to perform calculations 
in a two dimensional parameter space (lattice and spin 
disorder) , it is desirable to introduce the effect of temper- 
ature using more realistic models for the lattice and spin 
disorder. The method we have used can be straightfor- 
wardly extended in this direction as well as to more com- 
plex materials. Finally, we hope that this work will stim- 
ulate more experimental temperature-dependent studies 
of magnetization damping. 
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Appendix: A-k in Gilbert damping frequency 

The damping A has units of s,"^ so its numerical value 
should not depend on whether SI or Gaussian units are 
used. Conversion of the dimcnsionlcss parameter a to the 
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damping frequency A should follow the general relation 

A = 7 • a • magnetization density. (A.l) 

So converting from SI to Gaussian units should just re- 
quire converting the magnetization density, Mg (mag- 
netic moment per unit volume that we calculate as 
/LtB/atom) from SI to Gaussian (cgs) units. ^^'-^^ In SI, 
Ms is measured in units of A/m and the damping A is 
related to a by 

A = 7aAfs(SI). (A.2) 

In Gaussian units, the magnetization should be measured 
in units of magnetic moment cm"'^. If measured in Oer- 
steds or Gauss, then the magnetization is AirMg and A is 
related to a by 

A = 7 • a • AttMs (in Gaussian units). (A. 3) 



In most experiments, the magnetization is measured and 
reported in Gaussian units but Eq. A.2 is used instead of 
Eq. A. 3 which leads to a factor An missing in the damping 
frequency. 



For example, a recent measurement"^^ reports the mag- 
netization of bcc iron to be AirMg = 21.1 kG, corre- 
sponding to 1.68 X 10^ A/m in SI or 2.13 /ie/atom. 
Choosing g = 2 (for simplicity), we have 7 = 0.2203 
MHz-m/A= 17.588 MHz/Oe. We convert the reported 
a = 0.0019 to A and obtain the same frequency 7.1 x 10* 
Hz in SI using Eq. A.2 or in Gaussian units using Eq. A. 3. 
However, the frequency 57 ± 3 MHz reported in Ref. 31, 
which is obtained using Eq. A.2 combined with magne- 
tization in Gaussian units, has to be multiplied by Att to 
be consistent. 
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